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Massless neutrino exchange leads to a new long-range force between mat- 
ter. Recently, it was claimed both that the potential energy due to this 
interaction i) dominates the total energy of neutron stars and ii) that it 
is zero. We recalculate the energy of a neutrino propagating in a classical, 
uniform background of neutrons and find a negligible, but non-zero con- 
tribution to the total energy of neutron stars. We estimate the neutrino 
pair creation rate of a neutron star caused by a density gradient of the 
background neutrons but found it too small to be observable. 
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1 Introduction 



Recently Fischbach calculated the energy difference AE" between a neutrino 
immersed in a neutron star and in vacuum . He used perturbation theory to 
derive the potential energy VT^'^) of N neutrons due to the exchange of massless 
neutrinos, Ai? = The ratio of the contributions to /S.E from k and 

k + 2 body interactions found by him equals for k <^ N 



1 ^GpN\ 



{k + 2){k + l) V 



where is the Fermi constant and R the radius of the star. For a typical 
neutron star, GpN/R^ = O{10^^) and multi-body effects become dominating. 
In particular, |Aii^| exceeds the mass energy of the neutron star — an obvious 
contradiction to the observation of neutron stars. The resolution of this para- 
dox proposed in Ref. is to consider a massive neutrino. Then the neutrino 
can interact only with neutrons within its Yukawa radius and, if m,^ is 

sufficiently large, the dangerous many-body effects are exponentially damped. 
Hereby, a lower bound for the electron neutrino mass, m^^ ^0.4 eV, was 
derived 
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The method and results of Ref. [|T| provoked some criticism. Smirnov and 
Vissani ^ argued that a neutrino sea inside the neutron star |^ reduces the 
potential energy W^^'^ because of Pauli blocking. They stressed also that the 
behaviour of the potential energy W^^\ 

iy(2fc) = (_i)fc|H/(2fc)|^ (2) 



is unacceptable. In fact, both the increasing of iTy^*^-*! and its oscillatory be- 
haviour with 2k are clear signs for the breakdown of perturbation theory. 
Abada et al. Q recalculated the energy difference /\E taking into account 
non-perturbatively the interaction of the neutrino with a uniform neutron 
background. They obtained = if no neutrino sea is presented. This 
result would imply that neutrinos do not interact at all with the neutron 
background. In view of the results of Ref. and the possible implications 
for neutrino physics, we feel it appropriate to recalculate once again ^E. Ad- 
ditionally, we estimate the neutrino pair creation rate of a neutron star caused 
by a density gradient of the background neutrons. 

In our calculations we always assume that the neutron density n7v(x) respec- 
tively its gradient Vn7v(x) can be approximated locally by a constant value. 
This assumption looks unproblematic because of the macroscopic size of the 
neutron star, neglects however possible long-range effects due to massless neu- 
trinos. 



2 Neutrino energy density without neutrino sea 



The energy E oi a. neutrino interacting with a classical background current 
of neutrons is given by 

E = (P°) = ijd^x {0\ilj\x)dtip{x)\0)j . (3) 



Here, ^p is the field operator in the Furry picture obtained after second quan- 
tizing the solutions of the Dirac equation 



7^9'^ + im - V2iGFaNl^.J''PL\ H^) = 



(4) 



where qn = —0.5 and = (1 — 7^)/2 projects out the right-handed compo- 
nent of the neutrinos. 

Since we are interested in the infrared regime, we can assume in the following 
as static, = {nN,0). Moreover, the neutron number density un is also 



2 



nearly homogenous. Therefore, the solutions of Eq. ) can still be character- 
ized by the four-momentum = {E, p) of the neutrino. Then, as it is well 
known, the only change compared to the vacuum case is the modification of 
the dispersion relation of the neutrino [^|[ , 



E = {m^ + \p\^Y/^±V 



V 



'V2 



UN 



(5) 
(6) 



Here, the upper sign in Eq. (^) corresponds to neutrinos and the lower sign to 
antineutrinos. 



Equation can be rewritten as 

J x'\x 



(7) 



where Sf{x, x') is the Greens function of Eq. (^). Now we evaluate the energy 
difference 



= / d^x 



dp^ 

W) 



between a neutrino propagating in a neutron background |T(]| 

1 



Sf{p) 



ijp — V)'y^ — p ■ 7 — m + ie: 



Pl 



(9) 



and a neutrino propagating in the vacuum. 



sPip) 
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■/) — m + ie 



Pl. 



(10) 



Performing first the x and then the p integral results in 

Ai^ = -^/dpVV'^''(*-*')tr{/[SH/,0)-4°H/,0)]}^,^^. (11) 



Although each of the two contributions to the energy difference AE is UV- 
divergent, a finite final result for AE can be obtained combining the two 
fractions before integrating. Performing first the subtraction regularizes the 
integral, because then its leading term in p° vanishes, 



iV 



AE = -— I dp>V 

TT 



[(po - Vy -m? + is] [(p0)2 -m? + is] 



.(12) 
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The final integral can be done with the help of the residue theorem (cf. Fig. 1) 
and gives independently of the neutrino mass m 



AE = V^-19eV — (13) 
0.3 fm-3 ^ ' 



This result could be expected in virtue of Eq. (|^). Since AE is the energy 
difference of one neutrino, we obtain the total contribution AM of neutrinos 
to the self-energy of a neutron star with radius R ^ 10 km as 

4-7r 71 

AM = —R^n^AE ^ -42 kg ^. (14) 

3 ^ 3 X 10-22 fm^3 ^ ^ 



Obviously, AM is negligible for a realistic value (cf. next sections) of the 
neutrino density n^. 

Finally, we want to comment on the possible source of error in the results 
of Ref. [|l],3. In Ref. |l[], the 8th order term of perturbation theory W^^^ was 
used to estimate AE. Since J2k W^''^ is divergent and alternating, any result 
obtained in a finite order of perturbation theory is meaningless. By contrast, 
the authors of Ref. 0] used the same non-perturbative method as we did. 
However, they made - as apparent from their Eq. (8) - the limit x' \ x 
before integrating over x. Then, they argued that after regularization their 
d^p integral vanishes due to the antisymmetry of the integrand. However, this 
is not true because of the presence of the factor p'^e~^P \ 



3 Neutrino energy density with neutrino sea 



Loeb proposed first that neutrinosQ with energies below ~ 50 eV are bounded 
inside neutron stars, while antineutrinos are repelled 0. Qualitatively, this 
result follows directly from Eq. (^ considering a neutron star as potential 
wall with depth V and radius R, and simply assuming that all the levels with 
energy [0, V] are occupied. However, one should note that neutrino states with 
energy V < E < —niy do not have an exponentially damped wave function 
for r > R, and are resonances instead of true bound states. 

To account for a possible neutrino sea inside an neutron star, we can apply 
either the imaginary or the real-time formalism of finite-temperature field the- 
ory. We use the latter since it is easier in this formalism to separate the medium 

^ In Ref. |9|, it is argued that only Dirac neutrinos are trapped while Majorana 
neutrinos can leave freely the medium. 
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from the vacuum effects. Then the thermal neutrino propagator Sf{p,P) is 
given by 



p-7 



(po - vy 



ipi 



27ri5((/ - Vf - |pn/(/ 



Pt. 



(15) 



where /? = 1/T is the inverse temperature and f{p^) is the distribution func- 
tion of the neutrinos. To simphfy the notation, we have set m = 0. In the case 
of thermal equilibrium, f{p^) is the Fermi-Dirac distribution function 



9{p^ - V) e{-p^ - V) 



+ 



(16) 



Note that the edge between particles and antiparticles is shifted by V . Insert- 
ing the thermal propagator into Eq. (|^, performing first the x integration 
restricted to the volume V of the neutron star and then the po integral results 
in 



^Ep/V = I ^ [(IpI + V)U\P\ + V) + (-IpI + V)M\p\ - V)] (17) 

= / (^ [^./.(^.) - E,ME,)] (18) 
= {E,)p-{E,)p. (19) 

In the last step, we denoted the thermal average with the (anti-)particle dis- 
tribution functions by (. . Also in this case, the result obtained has a very 
plausible form. 

To obtain a numerical estimate for AEj^, we follow Ref. and assume a 
degenerated neutrino sea with rij, ~ 3 x 10~^^ fm~^ and Fermi momentum 
Pi? ~ 50 eV. Furthermore we set rip = and obtain 

AEs/V (v + -pf] ~ 5.5 x lO'^i ^ ^ . (20) 

V 4^ J fm^ 3 X 10-22 fm-3 ^ ' 



For a neutron star, the thermal contribution AEp is dominated by the contri- 
bution of the kinetic energy of the neutrinos and therefore positive. However, 
since the number density of neutrinos is much smaller than the number den- 
sity of neutrons, Uy/niq ~ 10-2^^ the thermal contribution AEp is irrelevant 
compared to AE. Finally, we want to remind that - as mentioned in the intro- 
duction - the results obtained in section 2 and 3 are only valid for an uniform 
neutron background, i. e. in the limit of an infinite neutron star. 



5 



4 Spontaneous neutrino pair creation 



Let us consider in more detail tlie analogy between a neutron star and a po- 
tential wall. We keep now again the neutrino mass rrii, finite and consider 
Dirac neutrinos. The limit for the electron neutrino mass is rriiy^ < 5 eV. 
Hence the potential V ^ —50 eV is overcritical, V < —rriu^, and there- 
fore able to produce t'e^'e-pairs at its interface. The quantity characterizing 
this process, m^/|V\^|, is not well-defined in the simple picture of the poten- 
tial wall. We assume instead that ni^{x) can be approximated by njv(x) = 
nAr(xo) + (x — xo)VnAr(x). Then we can treat locally the "electric" field 
E = —W — dtA, where A = (V,0), as uniform and derive the sponta- 
neous neutrino pair creation rate due to a density gradient of the neutron 
background. 



The action S{A) describing the vacuum with a potential A is given by |1L1 



In S{A) = Tiln [S^^SP] . (21) 



Here, Tr means the trace over Dirac indices and integration over the contin- 
uous variable x. Furthermore, Sp denotes the operator with matrix element 
{x\Sf\x') = Sf{x,x'). Adding the transposed version of its RHS to Eq. (PT]), 
inserting CC^^ = 1 and using C'-y^C~^ = 7^'*, Cj^C^^ = —•yjj^, we obtain 



2 In S{A) = Tr In S^Piif - 4 + - ie){f - ^ - + ie)PL 

1 



X 



+ is 



-Pl 



le 



-Tr In 
2 



(P _ Af -ml + ie + -a^.F^'" 



p2 _ ^2 _|_ 



(22) 



(23) 



This is 1/2 of the corresponding result for the pair production of fermions by 
an uniform electric field. Therefore, we can borrow the final QED result |l! 
and obtain for the pair creation probability per unit time and volume 



E 



2 00 



w 



E 



mrm 



87r3 — w 



exp 



\E\ 



(24) 



To obtain an order-of-magnitude estimate for w, we choose a density profile 
that models roughly a neutron star with a soft Reidl equation-of-state ]TB| : 

^ Since the absolute value of the potential V has no physical meaning, a more 
correct statement is that the potential difference AV = V{r < R) — V{r > R) is 
overcritical. 
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we assume as radius of the star R = 10 km, a homogenous core and a crust 
with thickness L = 2 km, in which the density decreases hnearly to zero. 



, uq for r < R — L 

nNir) = { (25) 

no{R-r)/L for R-L<r<R. 



Then, in a volume V ~ 2 ■ 10^^ cm^ of the neutron star exists a field \E\ = 
V/L ~ 5 • 10^^ eV^ in radial direction. If ml > \E\/{mT), the creation of 
a uu-paii by n field quanta is exponentially suppressed. But even if ^ 
(|E|/7r)i/2 ^ 4 . ^0-5 eV, the total luminosity C is only 

,11 , / 



£ = 2m,wV ~ 2.1 ■ lO^^erg/s ( ) (26) 



and has therefore practically no influence on the energy budget of the star. 
However, z/gi^e-pair creation by neutron stars is probably the only case in the 
present Universe that massive particles are spontaneously created and there- 
fore it is interesting in its own. 



5 Summary 



We have recalculated the energy of a neutrino in the uniform background of 
classical neutrons, both without and with a neutrino sea. In the first case, we 
found that the energy is changed by the small amount V ~ —20 eV, in agree- 
ment with the well-known result of Wolfenstein . Moreover, the influence 
of a possible neutrino sea inside the neutron star has an even smaller effect 
on the neutrino energy. Therefore, the contribution of neutrinos to the total 
energy of a neutron star seems to be negligible. However, to settle definitely 
the question if many-body effects become important in neutron stars it is nec- 
essary to calculate AE not only non-pertubatively but to take into account 



also the finite geometry of the star We have estimated the energy-loss of 



an neutron star due to spontaneous z/i/-pair production but found it too small 
to be observable. 
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After submission of this work, a preprint of As. Abada, O. Pene and J. 
Rodrigues-Quintero (|hep-ph /9 71 2266|) following the line of arguments of Ref. [§] 
appeared. Their main conclusion is that the border of a finite neutron star 
(modelled as a square wall) automatically generates the neutrino sea inside 
the star. 
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